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INTRODUCTION 


( 


^  A  confidence  ellipse,  soffletimes  called  an  EEP,  is  one  of  the  ways  in 
which  intelligence  data  concerning  an  emitter's  location  is  sunnarized.  Such 
.  an  ellipse  is  always  centered  at  the  estimated  emitter  location  a»'-strown-tn~^ 

}  rtgtire"t>  and  is  a  measure  of  the  possible  error  in  the  location  estimate. 

In  Figure  1 ,  X  is  the  estimated  location  of  the  emitter.  The  ellipse  is 
a  95t  confidence  ellipse,  tdiich  means  that  an  observer  may  be  95}  confident 
that  the  ellipse  contains  the  true  location  of  the  emitter.  Put  another  way, 
it  is  unlikely  that  the  estimated  location  is  in  fact  the  true  location.  The 
confidence  ellipse  defines  a  region  in  Uhlch  the  true  location  is  most  likely 
to  be  found.  Consequently,  the  size  of  the  ellipse  is  a  measure  of  how  good 
the  location  estimate  is. 

There  are  two  ways  of  specifying  the  shape  of  a  confidence  ellipse t 
geometrically  and  statistically.  The  first  method  is  that  used  in  Figure  1. 
The  ellipse  has  a  semi-major  axis  of  length  A  and  a  semi-minor  axis  of  length 
B,  and  is  rotated  a  degrees  from  true  north.  In  this  example  the  axis  lengths 
are  measured  in  kilometers.  The  second  method  is  not  as  easy  to  interpret 
geometrically,  but  la  mathematically  more  convenient.  Here  the  shape  and  size 
of  the  ellipse  la  defined  by  a  variance-covariance  matrix,  which  is  commonly 
denoted  by  the  Greek  letter  Z.-  Since  this  symbol  is  also  commonly  used  as  the 
mathematical  summation  operator,  the  variance-covariance  matrix  will  be 
denoted  by  the  letter  S.  S  is  generally  assumed  to  be  known  for  any  given 
sensor  system. 

Once  an  analysis  center  has  received  data  in  this  form,  a  step  called 
self -cor relation  is  commonly  performed.  The  goal  of  self-correlation  is  to 
'  determine  if  two  or  more  pieces  of  data  actually  concern  one  individual 
emitter.  It  is  clear  why  this  is  desirable. 

The  purpose  of  the  analysis  center's  data  base  is  to  store  information 
about  various  intelligence  objects.  If  two  pieces  of  data  come  in,  both  of 
which  concern  one  object,  such  as  an  emitter,  it  makes  sense  to  pool  the 
Information  in  some  way.  Increased  information  should  give  a  better  estimate 
of  the  emitter's  location.  Self-correlation  tests  whether  two  pieces  of  data 
should  be  pooled,  or  combined,  and  if  so  performs  the  combination.  Figure  2 
displays  a  simple  description  of  the  self-correlation  process.  Note  that  this 
is  not  intended  to  describe  any  particular  intelligence  system,  but  is  merely 
an  explanation  of  the  general  concept. 

A  new  location  estimate  and  confidence  ellipse  is  received  by  the 
analysis  system.  It  is  compared  to  each  ellipse  already  in  the  data  base.  If 
the  statistical  test  used  by  the  system  Indicates  that  the  new  data  refers  to 
the  same  emitter  as  an  ellipse  in  the  data  base,  the  new  ellipse  is  combined 
with  the  old  ellipse,  and  the  combined  data  is  stored.  If  the  statistical 
test  indicates  that  the  new  data  do  not  refer  to  an  emitter  in  the  data  base, 
the  new  data  is  entered  as  a  separate  entry. 

Figure  3  shows  two  ellipses  that  the  test  rejects  as  being  from  the  same 
emitter.  Figure  4  shows  two  that  the  test  accepts,  as  well  as  their  combined, 
or  resultant,  ellipse. 


^The  goal  of  this  paper  is  to  examine  the  geometric  implications  and 
properties  of  the  statistical  teat  for  deciding  whether  two  estimates  refer  to 
the  same  emitter,  and  of  the  method  of  combination  for  confidence  ellipses. 
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First,  understanding  these  properties  provides  an  intuitive  basis  for  looking 
at  self-correlation.  Will  two  ellipses  be  accepted  as  coming  from  the  same 
emitter?  How  will  they  combine?  Questions  such  as  these  may  frequently  be 
answered  simply  by  looking  at  the  two  ellipses.  Further,  the  test  and 
combination  method  are  based  on  a  number  of  assumptions  which  may  or  may  not 
be  met  in  actual  use.  If  one  of  these  assumptions  is  not  met,  then  the 
statistical  test  or  the  combination  algorithm  may  give  spurious  results.  The 
geometric  properties  described  in  this  report  provide  a  means  of  understanding 
how  robust  these  methods  are  with  respect  to  various  assumptions  —  that  is, 
how  large  an  effect  relaxing  an  assumption  will  have.  Lastly,  the  properties 
discussed  in  this  report  shed  light  on  the  kind  of  problems  vrtiich  will 
naturally  arise  >^lle  using  the  statistical  test  and  combination  method  for 
se 1 f -c or re la t ion . 

To  this  end,  the  following  points  will  be  discussed: 

Shape,  size,  and  orientation  of  the  resultant  ellipse. 

Location  of  the  resultant  ellipse. 

Geometrical  properties  of  combining  more  than  two  ellipses. 

Geometrical  properties  of  the  statistical  test  between  two  ellipses. 

Properties  of  the  statistical  test  used  on  a  sequence  of  ellipses. 


II 


CONCLUSION 


The  geometrical  properties  of  the  statistical  test  for  deciding  if  two 
estimates  come  from  the  same  emitter  and  of  the  combination  method  for 
confidence  ellipses  will  be  discussed  in  sections  III  and  IV.  However,  there 
are  certain  broad  results  %ihich  will  be  stated  here,  particularly  as  regards 
combining  or  testing  a  sequence  of  estimates. 

If  the  combination  method  is  used  to  combine  a  sequence  of  estimates,  the 
order  is  not  important;  the  final  resultant  ellipse  will  be  unchanged. 
However,  the  statistical  test  is  not  independent  of  order. 

Refer  back  to  Figure  2,  concerning  the  self-correlation  algorithm.  In 
Figures  5,  6,  7,  and  6  we  will  trace  a  series  of  ellipses  through 
self-correlation.  Note  that  the  ellipses  in  these  figures  (which  will  be 
called  A,  B,  and  C)  are  drawn  to  scale  —  if  the  text  says  that  the  test  will 
accept  two  ellipses  as  coning  from  the  same  emitter,  this  is  the  case  for  the 
two  ellipses  shown.  Where  two  ellipses  are  combined,  the  resultant  ellipse 
will  be  called  G,  H,  or  I. 

Case  1:  A,  B,  C. 

Initially,  there  are  no  ellipses  in  the  data  base.  Sensors  in  the  field 
send  ellipse  A  into  the  analysis  center,  where  it  goes  through  self- 
correlation.  Since  there  are  no  other  ellipses  to  compare  A  with,  it  is 
entered  into  the  data  base. 

Now,  ellipse  B  is  sent  in  from  the  field,  and  goes  through  self- 
correlation.  It  is  compared  with  A,  and  the  statistical  test  accepts  that  A 
and  B  refer  to  the  same  emitter  ~  that  is,  they  are  both  estimates  of  the 
location  of  a  single  emitter.  B  is  therefore  combined  with  A.  A,  B.  and 
their  combination  is  shown  in  Figure  5a. 

Since  A  and  B  were  combined,  there  is  only  one  ellipse  in  the  data  base: 
the  combination  of  A  and  B,  which  is  called  G.  G  is  shown  in  Figure  5b. 

A  third  ellipse  is  received  from  the  field,  ellipse  C.  It  is  compared 
with  G,  and  the  test  accepts  that  they  refer  to  the  same  emitter.  C,  G,  and 
their  combination  is  shown  in  Figure  6a.  The  combination,  ellipse  H,  replaces 
ellipse  G  in  the  data  base  (Figure  6b}.' 

Thus,  when  ellipses  A,  B,  and  C  are  received  and  tested  in  that  order, 
the  final  result  is  a  single  ellipse,  ellipse  H. 

Case  2:  A,  C,  B 

Now,  suppose  these  ellipses  were  received  in  a  different  order.  This 
could  easily  change  the  results. 

As  before,  A  is  the  first  ellipse  sent  in  from  the  field,  and  is  thus  the 
first  ellipse  in  the  data  base.  Next,  ellipse  C  is  received.  The  statistical 
test  rejects  that  A  and  C  refer  to  the  same  emitter,  and  so  no  combination  is 
done.  As  indicated  by  Figure  2,  ellipse  C  is  entered  into  the  data  base.  A 
and  C  are  shown  in  Figure  7a. 

Next,  ellipse  B  is  received  from  the  field,  and  is  compared  to  both  A  and 
C  (Figure  7b).  .The  statistical  test  accepts  that  A  and  B  refer  to  the  same 


emitter,  end  that  C  and  B  refer  to  the  same  emitter,  even  though  It  would  not 
accept  that  A  and  C  refer  to  the  same  emitter i  It  is  unclear  what  should  be 
done  in  such  a  case.  B  can  not  be  combined  with  both  A  and  C,  and  A  and  C 
cannot  be  combined  with  each  other. 

We  assume,  arbitrarily,  that  B  is  matched  with  ellipse  C  rather  than  A. 
Thus,  B  and  C  are  combined,  producing  ellipse  Z,  shown  in  Figure  8a. 

Finally,  ellipse  1  replaces  B  and  C  in  the  data  base  (Figure  8b). 

When  the  ellipses  are  received  and  tested  in  the  order  A,  C,  B,  the 
self -correlation  data  base  ultimately  contains  two  ellipse,  A  and  1.  Note 
here  that  if  B  had  been  combined  with  A  instead  of  C,  two  ellipses  still  would 
have  resulted.  Figure  8b  also  shows  ellipse  H,  the  result  of  Case  1 ,  for 
purposes  of  comparison. 

It  should  be  noted,  however,  that  the  exaople  above  is  based  only  on  one 
possible  version  of  the  self-correlation  algorithm.  Here,  an  ellipse  estimate 
is  only  compared  with  the  ellipses  in  the  data  base  when  it  is  first  received 
from  the  field.  A  more  sophisticated  self-correlation  process  might  also 
periodically  test  all  of  the  ellipses  in  the  data  base  with  each  other, 
combining  tdiere  appropriate.  If  such  a  process  had  been  used  in  the  example 
above,  then  the  results  of  the  two  cases  above  may  have  been  somewhat 
different.  The  point  to  be  made  here  is  that  the  order  in  which  the  ellipses 
are  tested  and  combined  will  still  make  a  difference  in  the  final  results. 

Thus,  while  the  order  is  not  important  to  the  combination  method  in  and 
of  itself,  the  order  of  the  ellipses  can  have  a  tremendous  impact  upon  the 
statistical  test,  and  thus  on  the  self -correlation  algorithm  as  a  whole.  This 
is  perhaps  the  most  Important  result  in  this  paper. 
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GEOMETRIC  PROPERTIES  OF  THE  COMBINATION  METHOD 


In  this  section,  a  number  of  properties  of  the  combination  method  are 
listed.  The  figures  are  at  the  end  of  the  section.  Mathematical  details  and 
proofs  are  contained  in  the  Mathematical  Appendix.  Note  that  all  of  these 
results  are  dependent  upon  the  ellipses  being  constructed  at  the  same  level  of 
confidence  (l.e.  95S,  90},  etc.). 


A.  The  resultant  estimate  of  target  location  is  dependent  upon  two 
original  point  estimates,  and  also  upon  the  shape,  size,  and 
orientation  of  the  original  ellipse  estimates. 

Comment:  Errors  in  ellipse  shape  and  size  can  affect  target 
location  estimates. 

B.  The  size,  shape,  and  orientation  of  the  resultant  ellipse  is 
dependent  only  on  the  size,  shape,  and  orientation  of  the  two 
original  ellipses,  and  not  on  their  relative  positions.  (See 
Figures  9,  10,  and  13). 

Comment:  This  implies  that  the  distance  between  two  estimates  of 
target  location  does  not  affect  the  algorithm's  'certainty*  of  the 
resultant  estimate  of  target  location.  This  assumes  that  we  are 
correct  in  accepting  that  the  two  estimates  are  from  the  same 
emitter. 

C.  If  the  original  point  estimates  are  equal,  the  resultant  point 
estimate  will  be  equal  to  them.  Stated  equivalently,  the  resultant 
confidence  ellipse  will  be  centered  at  the  same  point  as  the 
original  ellipses.  (See  Figure  11). 

D.  If  the  original  point  estimates  are  equal,  that  is,  the  original 
ellipses  are  centered  at  the  same  point,  then  the  resultant  ellipse 
will  be  wholly  contained  within  their  intersection.  (See  Figures  11 
and  12). 

Figure  13  demonstrates  both  point  B  and  point  D.  In  13a,  the 
resultant  ellipse  is  contained  within  the  intersection  of  the 
original  ellipses.  In  13h  the  positions  (point  estimates)  of  the 
original  ellipses  have  been  moved,  but  the  resultant  ellipse  has  the 
same  size,  shape,  and  orientation  as  in  13a. 

E.  When  the  original  point  estimates  are  not  equal,  the  resultant 
ellipse  will  be  approximately  positioned  upon  the  intersection  of 
the  original  ellipses.  (See  Figure  14). 

Comment:  This  is  not  the  case,  of  course,  when  the  ellipses  do  not 
intersect.  However,  the  statistical  test  will  not  accept  two 
ellipses  for  combination  if  they  do  not  intersect,  so  this 


qualification  is  irrelevant.  (See  Section  IV,  E). 

If  the  original  ellipses  are  tangent,  but  not  otherwise 
intersecting,  the  point  estimate  of  the  resultant  ellipse  will  be 
located  near  the  point  of  tangency.  Further,  if  the  sizes  of  the 
original  ellipses  are  increased  by  a  factor  of  the  square  root  of  2 
in  such  a  case,  the  intersection  of  the  expanded  ellipses  will 
contain  the  point  estimate  of  the  resultant  ellipse.  (See  Figure 
15,  16,  17,  and  18). 

Comment:  This  result  is  interesting  mainly  in  light  of  point  G, 
Immediately  below,  since  the  statistical  test  will  not  accept  two 
tangent  ellipses  for  combination.  (See  point  E  of  section  IV). 

If  the  two  ellipses  are  not  tangent,  the  resultant  point  estimate 
may  be  located  geometrically  as  follows:  decrease  (or  increase)  the 
size  of  both  original  ellipses  proportionally  until  they  are 
tangent.  Then  the  point  estimate  may  be  located  as  in  point  F 
above.  In  practice  the  size  of  the  ellipses  would  only  be 
decreased,  since  the  statistical  test  will  not  accept  two  ellipses 
which  do  not  intersect.  (See  point  E  of  section  IV.) 

When  more  than  two  ellipses  are  combined  the  order  of  combination 
does  not  change  the  result.  Thus  the  only  difference  order  of 
combination  can  make  la  when  the  statistical  teat  is  used  (to  accept 
or  reject  that  data -are  from  the  same  emitter).  (See  Figures  19  and 
2< 


More  examples  with  a  common  center.  Note  that  the  ellipse 
angle  ms  ahanged  slightly,  significantly  affecting  the  size 
of  the  resultant  ellipse. 


Figure  11 


Note  that  only  the  length  of  tne  axis  of 
one  of  the  ellipses  is  changed.  The  resultant 
ellipse  changes  only  slightly  from  ]2b  to  12c. 


Fiaure  12 
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In  15a  we  see  that  the  resultemt  locati/^n  estimate  is  close 
to  the  tangency  point.  Figure  15b  is  VT  larger  and  the 
location  estimate  is  inside  the  intersection  as  it  will 
always  be. 


Figure  15 


See  Figure  1$«  for  the  basic  idea.  The  difference  here 
is  that  the  original  ellipses  are  not  tangent  but  if 
they  are  shrunk  proportionately,  the  point  of  tangency 
is  close  to  the  location  estimate. 


Figure  16 


GEOMETRICAL  PROPERTIES  OF  THE  STATISTICAL  TEST 


The  following  geometrical  properties  are  proven  in  the  Mathematical 
Appendix.  Figures  are  at  the  end  of  the  section.  Note  that  these  results 
assume  that  the  statistical  test's  level  of  significance  and  the  ellipses' 
level  of  confidence  are  equal  (i.e.  if  the  test  is  a  951t  test,  then  all  of 
the  ellipses  must  be  95}  ellipses). 


A.  The  absolute  orientations  of  the  two  original  ellipses  have  no 
effect  on  the  statistical  test,  only  their  orientations  with  respect 
to  each  other. 

Figure  21a  shows  two  ellipses,  and  Figure  21b  shows  the  same  two 
ellipses  rotated  an  arbitrary  number  of  degrees.  The  test  result  is 
the  same  in  either  ease. 

B.  The  units  of  measure  have  no  effect  on  the  test.  One  axis  could  be 
measured  in  meters  and  the  other  in  kilometers. 

C.  Order  of  combination  can  matter  when  three  or  more  estimates  are 
involved. 

This  is  an  extremely  important  result.  An  excellent  example  is 
contained  in  the  Conclusion  (section  II).  Also  see  the  "worst  case" 
result,  below. 

0.  If  either  of  the  ellipses  in  the  test  contains  the  center  (point 
estimate)  of  the  other  ellipse,  the  test  will  definitely  accept. 

Figure  22  shows  several  cases  of  this  result.  The  statistical  test 
would  accept  in  each  case. 

E.  If  the  ellipses  do  not  intersect  at  all,  or  only  at  one  point 
(tangency),  then  the  test  will  definitely  reject. 

Figure  23  shows  several  cases  of  this  result.  The  statistical  test 
would  reject  in  each  case. 

F.  There  are  cases  which  do  not  fit  into  either  D  or  E  above.  For 
these  cases,  more  complex  geometric  constructions  are  necessary  to 
determine  whether  the  test  will  accept  or  reject. 

Figure  24  shows  several  cases  in  which  neither  point  D  nor  point  E 
can  determine  acceptance  or  rejection.  24a  and  24b  show  cases  in 
which  the  test  would  accept,  and  24c  and  24d  show  cases  in  which  the 
test  would  reject. 
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The  "Worst  Case"  Result  of  the  Impact  of  Testing  Order 


As  stated  In  point  C  above,  the  order  of  testing  the  ellipses  can  have  a 
significant  impact.  This  result  shows  the  worst  possible  errors  that  can 
occur  for  three  ellipses.  First,  however,  we  must  define  two  pieces  of 


notation. 

Assume  that  we  have  three  ellipses.  A,  B,  and  C.  In  order  to  test 
whether  or  not  two  of  these  ellipses  refer  to  the  same  sensor,  a  test 
statistic  is  calculated.  The  exact  formula  of  this  statistic  is  shown  in  the 
Mathematical  Appendix,  but  here  we  shall  denote  it  by  T().  For  example,  if  A 
and  B  are  being  compared,  the  value  of  the  test  statistic  would  be  T(A,B). 

In  order  to  test,  we  choose  a  value  k  tdiich  will  give  us  the  desired 
significance  level.  If  T(A,B)Ak,  then  we  accept  that  A  and  B  refer  to  a 
single  emitter,  otherwise  we  reject. 

Suppose  that  T(A,B)Sk  when  the  test  Is  performed,  and  so  the  two  ellipses 
may  be  combined.  We  will  denote  their  combination  by  ABB. 

Three  further  observations  must  be  made.  First  of  all,  if  T(x,y}-0  for  any 
two  ellipses  x  and  y,  then  x  and  y  must  be  located  at  exactly  the  same  point. 
Secondly,  the  value  of  T(x,y)  indicates  In  a  statistical  sense  the  distance 
between  the  locations  of  ellipses  x  and  y.  T(x,y}*k  is  the  largest  value  of 
T(x,y}  that  will  allow  the  statistical  test  to  accept  x  and  y  for 
combination.  Finally,  if  T(x,y)*2k,  it  is  clear  that  the  test  will  not 
accept  X  and  y  for  combination. 

We  will  show  two  possible  "worst"  results.  Case  1  and  Case  2.  Both  of 
these  cases  show  what  may  happen  when  A,  B,  and  C  are  tested  in  different 
orders.  These  results  are  dependent  upon  a  mathematical  result  which  is 
proved  in  the  mathematical  appendix,  vmich  says  that  no  matter  what  order  the 
ellipses  are  tested  in,  the  sum  of  the  test  statistics  is  constant. 

Case  1 : 

The  shape  and  location  of  A,  B,  and  C  can  be  chosen  so  that  the  following 
is  true: 


Order  of  Testing/Combination 

Teat 

Results  A,  B,  C 


1st  Test  T(A,B)  •  k 


2nd  Test 


T(ABB,C}  - 


A,  C,  B 


T(A,C)  •  2k 


T(ASC,  B) 


These  statements  state  that  when  the  ellipses  are  tested  in  the  order  A, 
B,  C,  the  test  will  accept  each  time.  The  end  result  is  the  ellipse 
ABBBC.  However,  if  testing  is  done  in  the  order  A,  C,  B,  the  test  will 
reject  ellipses  A  and  C  for  combination,  even  though  ABC  and  B  must  have 
the  same  location  (T(ABC}  ■0). 


The  shape  and  location  of  A,  B,  and  C  can  be  chosen  so  that  for  order  of 
testing: 


Test 

Results 

Order  of  Testinj 

A.  B.  C 

{/Combination 

A.  C.  B 

lat  Test 

T(A.B)  -  0 

T(A.C)  -  k 

2nd  Test 

TCABB.C)  -  2k 

T(AOC.  B)  -  k 

These  results  state  that  A  and  B  have  exactly  the  same  location;  and  that 
when  testing  is  done  in  the  order  A,  B,  C  the  ellipses  ABB  and  C  will  not 
be  accepted  for  combination.  When  testing  is  done  in  the  order  A,  C,  B, 
however,  all  three  ellipses  are  combined. 

There  is  no  mathematical  preference  for  either  order  in  the  cases 
discussed  above.  That  is.  there  is  no  reason.  Just  from  the  math,  why  one 
should  prefer  to  test  in  the  order  A.  B.  C  over  A.  C.  B.  and  vice  versa. 

The  only  basis  for  making  this  decision  is  knowledge  of  UhLch  ellipses 
actually  ^ould  be  combined — which  would  make  the  test  unnecessary.  However, 
it  la  unclear  how  often  cases  like  this  occur.  While  the  errors  themselves 
are  significant,  the  chance  of  making  them  may  not  be. 


V.  SIMULATION  RESULTS 
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In  the  previous  sections  we  have  discussed  the  "goodness"  of  the 
combination  method  and  statistical  test  in  mathematical  terms.  For  example, 
we  have  shown  that  the  results  of  the  combination  method  are  not  affected  by 
the  order  of  combination;  any  number  of  ellipses  may  be  combined  in  any 
order  and  the  same  resultant  ellipse  will  be  gained.  However,  we  have  also 
shown  that  order  of  combination  can  have  a  serious  effect  on  the  results  of 
the  statistical  test. 

In  this  section  one  aspect  of  the  "goodness"  of  these  methods  will  be 
examined.  More  specifically,  we  wish  to  know  if  the  resultant  ellipse 
sometimes  misses  the  true  location  of  the  emitter,  and  if  so  how  often.  We 
will  use  simulation  to  study  this  question.  Simulation  Involves  running  an 
experiment,  such  as  finding  out  if  the  resultant  ellipse  contains  the  true 
location,  a  large  number  of  times,  so  that  the  chance  of  it  happening  may  be 
estimated.  Thus,  a  computer  will  simulate  data  coming  in  from  sensors  in  the 
field,  and  generate  many  confidence  ellipse  and  point  estimates.  The 
combination  method  and  statistical  test  will  be  applied  to  these  estimates, 
and  the  results  compiled. 

Before  the  results  of  the  simulations  can  be  discussed,  however,  the 
possible  errors  and  the  simulation  procedures  must  be  described. 

The  simulations  use  a  simple  model,  containing  either  one  emitter  or  two 
emitters  positioned  close  together,  and  two  sets  of  sensors,  A  and  B.  Each 
emitter  is  some  radio  source  such  as  a  radar  unit.  By  set  of  sensors,  we 
mean  an  array  of  sensors  working  together  to  get  a  fix  on  an  emitter.  The 
computer  generates  location  data  from  each  set  of  sensors.  If  there  are  two 
emitters,  then  sensor  set  A  Is  sighting  one  emitter  and  sensor  set  B  is 
sighting  the  other.  If  there  Is  only  one  emitter,  then  both  are  sighting  it. 

Each  set  of  sensors  gathers  data  about  the  location  of  the  emitter  it  is 
sighting,  and  sends  the  data  into  the  analysis  center  in  the  form  of  a 
confidence  ellipse  and  point  estimate.  In  the  analysis  which  follows  it  is 
the  ellipse  and  point  estimates  which  are  simulated,  as  opposed  to  the  actual 
sensor  readings.  The  estimates  are  compared  using  the  statistical  test,  and 
if  the  test  accepts  that  there  is  only  one  emitter  the  ellipses  and  point 
estimates  may  be  combined. 

There  are  a  nuitber  of  errors  which  may  occur  during  this  process.  First 
of  all,  the  statistical  test  may  decide  incorrectly.  That  is,  it  may  decide 
that  there  are  two  emitters  when  there  is  only  one  (Type  I  error),  or  it  may 
decide  that  there  is  only  one  emitter  ^en  there  are  in  fact  two  (Type  II 
error).  As  the  labels  TsT>e  I  and  Type  II  error  imply,  these  errors  are 
common,  and  may  occur  whenever  a  statistical  teat  is  used. 

Another  sort  of  error  involves  the  combined  estimate.  Suppose  that  there 
is  only  one  emitter,  and  the  statistical  test  accepts,  deciding  correctly. 
Then  the  combination  method  may  correctly  be  used  to  combine  the  two  ellipse 
and  point  estimates.  There  is  a  possibility  that  the  combined  ellipse  will 
not  contain  the  true  location  of  the  emitter.  Since  the  ellipse  is  used  to 
show  the  area  in  which  the  emitter  is  likely  to  be  found,  this  error  is  quite 
serious. 

A  more  complex  situation  results  if  there  are  two  emitters  but  the  test 
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decides  that  there  is  only  one  emitter.  In  this  case  one  would  hope  that  the 
combined  ellipse  would  contain  both  emitters,  in  order  to  avoid  compounding 
the  error. 

Thus,  we  are  interested  in  the  following  cases: 

A.  How  often  the  statistical  test  accepts  that  there  la  one  emitter. 

B.  If  It  accepts,  how  often  the  combined  ellipse  contains  the 
emltter(s). 

It  Is  simple  to  determine  mathematically  how  often  the  statistical  test 
should  accept  that  there  Is  one  emitter,  but  how  often  the  latter  case  occurs 
is  more  difficult.  Simulation  la  useful  here.  By  using  the  computer  to 
generate  the  data,  the  statistical  test  and  combination  method  may  be  applied 
many  times.  Each  time  one  of  the  errors  discussed  above  Is  made.  It  is 
recorded.  Thus,  the  probability  of  making  these  errors  can  be  estimated. 

The  first  step  Is  choosing  parameters  for  the  data  gathered  by  sensor 
sets  A  and  B.  The  program  assumes  that  there  are  two  emitters,  and  that  the 
data  gathered  by  each  set  follows  a  bivariate  normal  distribution  about  the 
true  location  of  the  emitter  being  sighted  by  the  sensor  set.  The 
variance-covariance  matrices  of  the  data  are  specified,  as  are  the  locations 
of  the  emitters.  The  case  in  which  there  Is  a  single  emitter  is  specified  by 
making  the  two  emitter  locations  equal.  Also  specified  are  the  number  of 
observations  each  set  of  sensors  will  use  to  calculate  its  confidence 
ellipses. 

In  the  second  step,  data  sets  of  the  specified  sample  size  are  generated, 
and  the  statistical  test  is  performed.  If  the  test  accepts  that  there  is 
only  one  emitter,  then  the  ellipses  are  combined,  and  it  is  determined 
whether  the  resultant  ellipse  contains  either  or  both  of  the  true  emitter 
locations.  In  the  results  discussed  later  in  this  section,  this  step  was 
repeated  500  times. 

Finally,  the  results  are  totalled  and  compared  to  the  total  number  of 
simulations  (here,  500).  In  this  way  probabilities  may  be  estimated.  For 
example.  If  the  test  accepts  475  times  out  of  500,  there  is  about  a  95% 
chance  that  the  test  will  accept  in  this  case. 

The  compiled  results  are  listed  In  Tables  1-4.  Consider  Table  1.  The 
parameters  for  the  runs  contained  in  Table  1  are  listed  in  the  heading. 

Sensor  set  A  generated  5  observations  about  the  true  location  of 
Emitter  1  at  (0,0),  with  covariance  matrix  for  each  simulation.  Also, 

Sensor  set  B  generated  5  observations  about  the  true  location  of  Emitter  2 
with  covariance  matrix  S^. 

The  Results 

The  first  column  of  the  table  show  the  location  of  Emitter  2.  When 
Emitter  1  and  Emitter  2  are  at  the  same  location  (0,0),  this  Is  equivalent  to 
there  being  only  one  emitter. 

The  second  column  shows  how  often  the  statistical  test  accepted  that 
there  was  only  one  emitter.  In  Table  1  this  ranges  from  97. 4S  of  the  time  to 
54.4}  of  the  time,  which  are  reasonable  values  given  the  parameters  listed  at 
the  top  of  the  table. 

Columns  three  and  four  show  how  often  the  combined  ellipse  does  not 
contain  the  true  location  of  Emitter  1,  Emitter  2,  and  both.  Colunn  three 
expresses  this  as  a  percentage  of  the  total  simulations,  and  column  four  as  a 


percentage  of  the  number  of  times  the  statistical  test  accepted.  For 
example'}  when  Emitter  1  and  Emitter  2  are  the  same,  in  4.6}  of  the  500 
simulations  the  combined  ellipse  did  not  contain  the  true  location  of  the 
emitter.  This  is  4.7}  of  the  total  number  of  times  the  statistical  test 
accepted  that  there  was  only  one  emitter.  (Remember  that  ellipses  are  only 
combined  if  the  test  accepts). 

Column  five  is  a  key  showing  what  is  in  each  block  of  percentages.  For 
example,  tdien  the  location  of  Emitter  2  is  (.5,  .5),  Emitter  1  is  not  in  the 
combined  ellipse  4}  of  the  time.  Emitter  2  is  not  in  the  ellipse  4}  of  the 
time,  and  neither  one  is  in  the  ellipse  2}  of  the  time.  Note  that  there  is 
double  counting  here;  for  exangile,  the  2}  of  the  time  that  neither  one  is  in 
the  ellipse  is  also  counted  in  the  4}  of  the  time  that  Emitter  1  is  not  in 
the  ellipse. 

Tables  2,  3  and  4  are  formatted  identically. 

Discussion  of  the  Results 

The  statistical  test  used  in  the  simulations  was  a  95}  test.  This  means 
that  when  there  really  is  only  one  emitter,  the  chance  of  saying  that  there 
are  two  is  limited  to  5}  (100}'95}).  As  stated  above,  this  is  called  Type  I 
error.  Indeed,  Tables  1 >4  Indicate  that  when  there  is  only  one  emitter,  the 
test  accepts  approximately  95}  of  the  time.  This  is  usual  and  reasonable. 
However,  there  is  an  additional  error:  sometimes  when  the  test  accepts,  the 
combined  ellipse  does  not  contain  the  true  emitter  location.  The  tables 
indicate  that  this  occurs  approximately  5}  of  the  time  as  well.  Since  this 
error  can  only  happen  when  the  test  has  accepted,  this  means  that  the  results 
are  in  error  about  10}  of  the  tlme—5}  for  this  and  5}  incorrect  rejections. 

This  may  be  a  significant  problem. 

Xt  is  also  enlightening  to  examine  what  happens  when  there  are  really  two 
emitters.  Since  there  are  two  emitters,  it  is  an  error  to  accept  there  is 
only  one  (Type  II  error).  It  is  better  to  reject,  and  deal  with  each  emitter 
separately.  However,  the  test  will  sometimes  accept  incorrectly,  and  combine 
the  confidence  ellipses.  From  a  practical  standpoint,  it  is  desirable  that 
the  combined  ellipse  contains  both  emitter  locations  in  this  case.  If  one  of 
the  locations  falls  outside  the  combined  ellipse,  the  analysis  center 
essentially  has  no  idea  of  its  whereabouts.  If  both  locations  are  in  the 
ellipse,  at  least  both  are  in  a  region  where  the  analysis  center  expects  to 
find  an  emitter. 

Examine  the  case  in  Table  2  where  Emitter  2  is  located  at  (1.5, 1.5). 

Here  there  is  approximately  a  55}  chance  of  committing  Type  II  error,  since 
that  is  how  often  the  test  accepted.  Note  that  in  almost  every  case  in  which 
the  test  accepted,  one  or  the  other  emitter  was  not  contained  in  the  combined 
ellipse  (88}  of  the  time).  About  half  of  the  combined  ellipses  contained 
Emitter  1,  and  about  half  contained  Emitter  2. 

Similar  results  to  these  may  be  found  in  each  of  the  tables. 

It  is  uncertain  how  often  this  problem  will  arise,  since  it  is  dependent 
upon  two  emitters  being  located  very  close  together.  "Close"  is  with  respect 
to  the  possible  spread  in  the  fixes  on  each  emitter's  location.  It  may  be 
that  signal  parametrics  might  help  in  distinguishing  between  them  as  well. 
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ELLIPSE  TEST  AND  COMBINATION  SIMUUTION  RESULTS 
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Simulations  •  500 
Sample  sizes  -  (5t  5) 

Using  TRUE  covariance  matrices 
Location  of  Emitter  1  ■  (0,0) 


V(’?  iS)-  Vf’o  io) 


Location 

of 

Emitter  2 

Accepts 

i  of  Total 

Accepts 

5  of  Total 

not  in 

5  of  Accepts 

Emitter  Location 
not  in 

Combined  Ellipse 

A. 65 

A.75 

1  Not  In 

(  0,  0) 

97. A* 

A. 65 

A. 75 

2  Not  In 

A. 65 

A.75 

Both  Not  In 

A.05 

A. 35 

1  Not  In 

{  .5.  .5) 

93. 0$ 

A. 05 

A.35 

2  Not  In 

2.05 

2.25 

Both  Not  In 

7.65 

8.35 

1  Not  In 

(  1.  1) 

92.05 

7.65 

8.35 

2  Not  In 

1.85 

1.95 

Both  Not  In 

11.85 

13.95 

1  Not  In 

(1.5.  1.5) 

8A.85  1 

13.65 

16.05 

2  Not  In 

\«e«/f  ^  9 10  f  ^ 

1.85 

2;15 

Both  Not  In 

1A.65 

19. A5 

1  Not  In 

(  2,  2) 

75.25 

19.25 

25.55 

2  Not  In 

0.85 

1.15 

Both  Not  In 

20.25 

29.75 

1  Not  In 

(2.5,  2.5) 

68.05 

25.05 

36.85 

.  2  Not  In 

3.25 

3.25 

Both  Not  In 

22.  A5 

Al.25 

1  Not  In 

(  3.  3) 

5A.A5 

27.25 

50.05 

2  Not  In 

3.85 

7.05 

Both  Not  In 

> 


• 

TABLE  2: 

4 

ELLIPSE  TEST  AND  COMBINATION  SIMULATION  RESULTS 

Sioulatloni 

1  -  500 

V^'o  10^ 

Saople  sizes  ■  (20,  20)  10^' 

Using  TRUE  covariance  matrices  ^ 

Location  of  Emitter  1  -  (0,0) 

Location 

of 

Accepts 

Accepts  not  in 

Emitter  Location 
not  in 

Emitter  2 

%  of  Total 

%  of  Total  %  of  Accepts 

Combined  Ellipse 

(1.5,  1.5) 


(  2.  2) 


(2.5,  2.5) 


1 

Not 

In 

2 

Not 

In 

Both  ! 

Not 

In 

1 

Not 

In 

2 

Not 

In 

Both  ] 

Not 

In 

1 

Not 

In 

2 

Not 

In 

Both  1 

Hot 

In 

1 

Not 

In 

2 

Not 

In 

Both  1 

Hot 

In 

1 

Not 

In 

2 

Not 

In 

Both  1 

Hot 

In 

1 

Not 

In 

2 

Not 

In 

Both  I 

Hot 

In 

1 

Not 

In 

2 

Not 

In 

Both  Not 

In 

K'' 
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Simulationi 
Sample  siz< 
Using  TRUE 
Location  ol 

Location 

of 

Emitter  2 

1  -  500 
ss  -  (5,  5) 
covariance  mati 
r  Emitter  1  -  (( 

Accepts 

%  of  Total 

•ices  *  u  c 

).0) 

Accepts  not  in 

K  of  Total  K  of  Accepts 

^b"^  0  10^ 

Emitter  Location 
not  in 

Combined  Ellipse 

i 

i 

•-VA 

^:.v 

4.6K 

4.8K 

1  Not  In 

(  0,  0) 

96. OK 

4.6K 

4.8K 

2  Not  In 

1  1 

i|.6K 

4.8K 

Both  Not  In 

9.4K 

10. IK 

1  Not  In 

(  .5.  .5) 

93. 4K 

7.8K 

8.4K 

2  Not  In 

*  *  ‘  *.• 

•  '  1 

2.-4K 

2.6K 

Both  Not  In 

19. 6K 

21 .8K 

1  Not  In 

Hss 

(  1.  1) 

89. 8K 

21 .8K 

24. 3K 

2  Not  In 

•  5.4K 

6.0K 

Both  Not  In 

36. 4K 

46. 5K 

1  Not  In 

(1.5.  1.5) 

78. 2J 

40.  OK 

5l;2K 

2  Not  In 

'i  • 

• 

13. 2K 

16.9K 

Both  Not  In 

;-.w 

47. 8K 

75. 9K 

1  Not  In 

(  2.  2) 

63.  OK 

50;  OK 

79. 4K 

2  Not  In 

34.8K 

34. 8K 

Both  Not  In 

m 

43. 2K 

93.1K 

1  Not  In 

'  . 

(2.5,  2.5) 

46. 4K 

42;4K 

91. 4K 

2  Not  In 

39.2K 

84.5K 

Both  Not  In 

28. 6K 

95. 9K 

1  Not  In 

MM 

(  3.  3) 

29. 8K 

29. 6K 

99.3K 

2  Not  In 

28. 4K 

95.3K 

Both  Not  In 

•:  I 

TABLE  4: 

ELLIPSE  TEST  AND  COMBINATION  SIMUUTION  RESULTS 
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Simulations  •  500 
Sample  sizes  •  (20,  20) 

Using  TRUE  covariance  matrices 
Location  of  Emitter  1  -  (0,0) 


Location 

of 

Emitter  2 

Accepts 

%  of  Total 

Accepts 

5  of  Total 

not  in 

5  of  Accepts 

Emitter  Location 
not  in 

Combined  Ellipse 

5.65 

5.95 

1  Not  In 

(  0,  0) 

94.45 

5.65 

5.95 

2  Not  In 

5.65 

5.95 

Both  Not  In 

23.45 

25.85 

1  Not  In 

(  .5,  .5) 

90.65 

20:85 

23;05 

2  Not  In 

5.85 

6.45 

Both  Not  In 

45.85 

76.15 

1  Not  In 

(  1,  1) 

60.25 

44.85 

74.45 

2  Not  In 

•  30.45 

50.55 

Both  Not  In 

31.85 

98.85 

1  Not  In 

(1.5,  1.5) 

32.25 

31.85 

98.85 

2  Not  In 

• 

30;  85 

95.75 

Both  Not  In 

9.65 

100.05 

1  Not  In 

(  2.  2) 

9.65 

9.65 

100.05 

2  Not  In 

9.65 

100.05 

Both  Not  In 

1.25 

100.05 

1  Not  In 

(2.5,  2.5) 

1.25 

1.25 

100.05 

2  Not  In 

i;25 

100.05 

Both  Not  In 

0.05 

— 

1  Not  In 

(  3.  3) 

0.05 

0.05 

— 

2  Not  In 

0.05 

— 

Both  Not  In 

MATHEMATICAL  APPENDIX 


DEFIHITIOWSt 

X  denotes  the  coordinates  of  an  arbitrary  point. 

X  denotes  the  coordinates  of  the  center  of  the  ith  ellipse  (l.e.  the  ith 

estimate  of  target  location.) 

denotes  the  covariance  matrix  of  the  ith  ellipse. 

denotes  the  set  of  points  in  the  1th  ellipse*{X| (X'X^}^S^^ (X-X^)$k} 

where  k  denotes  the  confidence  cut-off  value  being  used. 

0^  denotes  an  ellipse  similar  to  E^.  D^-{X|(X“X^)'^S“^ (X-X^)Sk/2) 

f  denotes  an  operator  on  the  X^,  S^,  and  E^  defined  so  that 

X^fXj-(S^«Sj ) X^+SJ^  Xj ) 

Ejf Ej-{X| (X-(X^eXj ) )“^ {X-(X^®Xj ) Xk) 

Note  that  •  is  a  mathematical  representation  of  a  method  of 
combining  ellipses. 

T(E^,Ej)  denotes  the  value  of  the  acceptance  test  statistic. 

T(E^,Ej)  - 


PROPERTY  1 ;  •  is  commutative. 


Proof.  SjeSj-(S"'  )'■*  -(SJ^  ♦S”’  )”'* -Sj®S. 
Cj.Xj.(Sj.Sj)(Sj'Xj.Sj'Xj).(Sj#SiKSj’Xj.SjTi,).Xj.Xj. 


PROPERTY  2;  •  is  associative.  (Together  with  Property  7  this  implies 
that  order  of  computation  does  not  matter.)  f  generalizes  to  combining  n 
ellipses  in  the  following  natural  way. 


i*t®i  usri  ' 


Proof.  Associativity  follows  directly  from  these  generalized  definitions, 
which  will  be  validated  by  induction.  When  n>2,  the  preceding  definitions  are 
the  original  definitions  for  t.  Now,  assume  they  are  valid  for  n*k.  It 
must  now  be  shown  that  they  are  valid  when  we  combine  the  k’«'1st  ellipse  with 
the  previous  combination  of  k  ellipses. 


Thus  it  is  true  for  covariance  matrices.  For  target  location  we  have 


llii=‘i-<i*!f’‘i>»=s.M 


.«,jfs,).v,)«,8?Si>-’<i!^i)*s-l,v,) 

■<5§lsi ' '  ‘jS?®i  >*'  «i8^i  >  <i5!f®r  *1  > 

k+1  k*1 

PROPERTY  3:  The  shape  of  the  resultant  ellipse  depends  only  on  the  shape  of 
the  original  ellipses  (Note  that  the  shape  of  any  ellipse  E  depends  only  on 
the  shape  of  its  covariance  matrix  S). 

Proof.  The  formula  for  S^fS^  depends  only  on  and  and  not  on 
or  Xj. 

PROPERTY  »;  If  X^-Xj  then  X^«Xj-X^. 

Proof.  XjiX  .(S,»S,)(S7'x,*ST'x,) 

-<S7'.sj’)-'(s7'Xj.Sj'Xj) 

■«r  *8]' ‘Sj'iXi 

■*i  :  • 

PROPERTY  5;  If  then  lies  entirely  inside  the  intersection 

of  E^  and  Ey 


Proof.  X  c  E^eEj  ->  {X-X^»Xj)^(SjtSj)“^(X-Xj«Xj: 

->  (X-Xj)^(S"NsJ^)(X-X^)gk 


i  _  r  J  ''  i' 

->  (X-Xj)^S^HX-Xj)Sk 
->  X  e  E^ 

Thus  Is  contained  in  E^.  Containment  in  E^  is  similar. 

PROPERTY  6t  If  X^-Xj  then  E^tE^  contains  the  intersection  of  and  Dj. 

Note  that  the  D  ellipses  are  proportionally  smaller  (.707)  versions  of  the 
E  ellipses. 

Proof.  X  c  Dj  ->  (X-X^)'^Sj^(X-Xj)Sk/2 

->  (X-X^iXj)^S"^(X-X^«Xj)Sk/2 
Similarly  X  c  ->  (X-XjiXj)'^Sj^  (X-X^tX^)Sk/2 

If  X  is  in  the  intersection  of  and  0^  then  X  c  0^  and  X  c  Dy  Since  both 
inequalities  hold  in  this  oase  they  may  be  added  to  get  that  X  satisfies 

(X-X^fXj )^(S'^ ♦Sj’ ) (X-XjfXj )Sk/2^k/2 
->  (X-X^iXj )^(S^fSj )"’ (X-X^*Xj )gk 
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->  X  c  E^Mj. 

PROPERTY  7t  is  »  point  Of  tangency  between  (X-X^)^S“^ (X-Xj)»Oj 

end  (X-Xj)TsJ^(X-Xj)-Oj  where 

o^-(X^eXj-X^)^S*^  CX^«X^-Xj, ) 

Proof.  The  gradients  of  the  two  curves  respectively  are 
GRAD j -2S*^ ( X-X^ )  GRADj -2S ( X-Xj ) 

At  XjgXj ,  GRAD^-2Sj^ ( (S^GSj ) ( X^+sj’ ) 

-2S^'  (  (S“^  ♦Sj^  (S“^  Xj^Sj^  Xj  )“Xj ) 

■2S-'  (S-'.S-' )-' ((S7'Xj.Sj'Xj)-(s;'.Sj’ )Xj) 

-2S“1 (ST^ ♦sj^ )'^  sj^ 

-SCSjS^^S  .S,Sj'Sj)-’(Xj-X^) 

-2(Sj+Si) 

Slallarly  at  X^gX^  GRADj-2(S^+Sj )"! (X^-Xj)— GRAD^.  Hence  tangency. 

PROPERTY  6;  X^gXj  la  the  point  where  (X-X^)‘^S“^  (X-X^)»{X-Xj)TsjUx-Xj) 

Is  a  minimum. 

*  Proof.  The  derivative  of  this  expression  using  the  notation  of  the 

previous  proof  Is  GRAD^’»’GRAOj .  The  previous  proof  shows  this  sum  Is  0. 

Note  that  since  the  Joint  likelihood  Is  a  function  of  the  quadratic  form 
minimized  In  PROPERTY  8,  It  Is  Implied  that  •  yields  a  maximum  likelihood 
estimate  of  the  associated  normal. 

PROPERTY  9;  Given  a  such  that  (X-Xj)Ts“1 (X-X^)»o  and  (X-Xj)Tsj’(X-Xj)-o 

are  tangent,  then  X-X.fX,  satisfies  (X-X. (X-X,  )S2o  and 
(X-Xj)Tsjl(X-Xj)S2a.  ^  ^  ill 

Proof.  By  PROPERTY  8  the  value  of  (X-X^)TsJ^ CX-Xj)^(X-Xj)TsJ^ (X-Xj)  can 
not  be  smaller  at  the  point  of  tangency  than  it  is  at  X^tX^  where  the 
value  la  a’*^a  or  2a. 

DEFINITION:  For  any  change  of  coordinates  of  the  form 

•(X)-AX>Y 

where  A  Is  any  Invertible  matrix  and  Y  Is  a  constant  vector,  define 
•(S}"ASA^  for  covariance  matrices  S. 

^  PROPERTY  10:  Change  of  coordinates  have  the  following  properties: 

(1)  ♦“^(X)-a”'x-a“^Y  and  ♦"' (S)-a"^S(a“' )^. 


(2)  la  'preserved*  under  this  change  of  coordinates,  l.e. 
lX|(X-Xj)Ts'^(X-X^)Sk}-{Xl(*(X)-*(X^))T*(S^)"^*(X)-*(Xj))ik} 

(3)  ♦(SjtSj)-*(S^)»*(Sj). 

(i»)  ♦(X^tXj)-*(X^)f*(Xj)! 

(5)  The  test  to  accept/reject  combination  of  and  Ej  Is  Independent 
of  the  coordinate  system  It  Is  calculated  in. 

Proof.  (1)  It  suffices  to  show  that 
For  X  note  that  (X))-*(il“^X-A"‘*y)-A(A"^X-A“^y)+y-X-Y+y-X. 

For  S  note  that  ♦(♦“' (S))-*(A‘^S(A"b^)-A(A'’s(A"b^)A^-S. 

(2)  (♦(X)-t(X^))T#(Sj)"^(*(X)-*(Xj)) 

-( ( AX*Y)-(AXj*Y)  )^(AS^a‘^)“^  (  <AX+Y)-(AXj^y)  ) 
-(A(X-X.))‘r(A^)“^ST^A“^(A(X-Xj)) 

-  ( X-X^  )^a'^  ( A^  )  S“ '  A“U  ( X-X  j  ) 

-(X-X^)Ts”^(X-X^). 

(3)  *(3^18  )-A(S7^+STb“U^ 

-((ATri'sTUsTb-U-b-^ 

.((AT)-is:V^:(AT)-isTU-b-‘' 

m  _e  tm  _4  ^4 


(4)  ♦(Xj«Xj)-A(X^eXj)>Y 


.((ASjA^)"^*(AS.A^)”^)"^ 


.A(Sj.Sj)(S-’Xj<Xj).I 

-A (3. fS. ) (3?^ Xi*ST’ X . )+A(S  fS . ) (3  f S. A"^  Y 
-AC3^«Sj)AT(AT)-^(s7^X,.Sj^Xj:(3"t3‘bA“'Y) 
-♦<3^»Sj)(a'^)"^  (37^"“’  (AX^*Y)+SjV''  (AXj+Y)) 
-♦(3^)f*(3j)((A3^A^)"^*(Xj)+(A3jA^)"’*(Xj)) 
-♦(3j  )  (♦(3^  )”^  ♦(X^  ) ) 

-♦(X^)f#<Xj). 


The  fact  that  properties  of  Interest  are  preserved  under  transformations 
tf  ths  type  specified  by  •  above  shall  be  used  to  simplify  the  analysis  in 
the  proofs  that  follow.  For  more  details  on  the  relationships  between  these 
transformations  and  symmetric  quadratic  forms  (such  as  ellipses)  and  the 
underlying  geometry  see  Metric  Affine  Geometry  by  Snapper  and  Troyer.  Some 
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readers  may  have  encountered  this  material  In  connection  with  the  study  of 
the  mathematics  of  Minkowski  space  and  it  is  the  pages  Immediately  proceeding 
this  example  that  are  of  the  most  Interest  here.  In  general  a  quadratic 
form  may  be  transformed  into  a  diagonal  matrix  of  1's,  O's  and  -1's.  Since 
the  quadratic  forms  of  interest  are  ellipses  only  1's  are  possible. 

For  the  purposes  of  proofs  one  of  the  center  points  may  be  sent  to  the 
origin  and  one  of  the  covariance  matrices  made  into  the  identity  matrix  by  use 
of  these  transformations.  The  one  qualification  is  that  the  results  that  are 
demonstrated  must  be  transformed  back  into  the  original  apace  if  shape  is  used 
in  the  result. 

PROPERTY  1 1 ;  If  the  center  of  one  ellipse  is  in  the  second  ellipse  then 
the  acceptance  test  will  accept.  On  the  other  hand,  if  the  ellipses  do  not 
intersect  at  more  than  one  point  then  the  test  will  reject. 


Proof.  For  the  first  part  of  the  proof  the  apace  will  be  transformed  so 
that 

x,.(o.o)  Xj.(..o)  s,.(;  0)  Sj.(;  «). 

Note  that  has-been  transformed  into  a  circle. 

First  it  must  be  shown  that  X^c  Ej  ->  (X2-X, )^(S^+S2)"! {X^-X2)Sk 
{X^-X2)‘^(S^+S2)"^  (X^-X2)-a2(d'H  )/((b+1  Xd^l  )~c^) 

Sa2d/(bd-c2) 

•  (X^-X2)‘rS2^(X^-X2) 

Sk  (since  X^  c  E2.) 

Consequently,  the  hypothesis  test  will  accept. 

For  the  second  part  of  the  proof,  the  space  will  be  transformed  as  follows, 
without  loss  of  generality. 


Further,  the  two  ellipses  are  assumed  to  be  tangent  (intersecting  only  at  one 
point),  tAere  the  point  of  tangency  P  is  as  follows: 

where  k  is  the  cutoff  value  defined  above. 

The  normal  vector  to  an  ellipse  at  X  is  of  the  form  S~^(X^-X}. 

The  normal  vector  to  at  point  P  is  S“^(X^-P)  •  I(O-P)  -  -P,  which  is 
proportional  to  the  vector  (^}.  Since  and  E2  are  tangent  at  P,  the 

normal  to  E2  is  also  proportional  to  (^},  and  is  thus  perpendicular  to  the 
vector  (0).  Since  the  normal  to  E2at  P  is  S2^(X2-P}|  we  know  that 
(0,  1 )SJ’ (X2-P)-0. 

Thus  (b,  c)(X2-P)*0,  and  Xj-P  is  proportional  to  the  vector  C"^).  This 
ioplies  that  X2«P’»(~^)t  for  some  t. 

We  find  the  value  of  t  as  follows,  letting  D  ■  Det  S2: 


-  ' 


.V 'A'.V. I*.  A  A ;■'. L»< L m ■  ^ '■•n H .' ■ '.' »j.  ■  j.' » ). ■^' '». j f.i  f> m  v» 'j  V >y. .'j.v* ■ 
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*■  •** 


(X2-P)^Sj^(X2-P)-k 
->  t2(-g)T(J  J)(-C)-k 

->  t2(-D,  0)("®)-k 

Solving  for  t  results  In  t-ivTcTTcOT.  The  negative  root  is  the  one  with 
the  correct  geometry  for  exterior  tangency. 

The  teat  will  only  accept  if  {X,-X2)^(I*S2)“Mx^-X2)Sk.  Substituting  in 
known  Information  results  in  (-X2)^(I*S2)“^ ("X2)Sk. 

Now  -/K  may  be  factored  out  of  X2  both  times,  allowing  the  k  on  the  right  to 
be  cancelled.  Thus  the  test  will  accept  if 

Simplifying  the  above  inequality,  we  get 

^  2  , 

<->  “JC'®®))S(a^D)(c^D)-b2 

<->  D(a+D)V57e(-D-cD)-v'57c(-ca-cD>b2)+{l/c)(-c(-D-cD)+b2D)S(a>D)(e+D)-b2 

(In  making  the  last  simplification,  remember  that  D  -  Det  S2-  ac-b  ). 

<->  D  ( a-^D ) +2/0/0  ( D+cD)  ♦  ( 1  /c )  C  c  ( c+ 1 )  ♦b^  ]DSac+ ( a+c )  D+D^-b^ 

<•>  aD^2/D/c(c+1 )D*(1/c)(c^+b^)DS(a+c)D 
<•>  a+2/D/c(c+1 )*(l/c)(c2+b2)Sa+c 

Finally,  the  test  will  accept  if  the  ellipses  have  an  exterior  tangency  and 
(b2/c)+2^D7c(c*1)S0. 

Since  c  and  D  are  always  greater  than  0,  this  inequality  is  never  satisfied. 
Thus,  if  the  ellipses  have  an  exterior  tangency,  the  test  will  reject.  From 
this  it  follows  that  if  the  ellipses  have  no  intersection,  the  test  will  also 
reject. 

The  next  result  is  not  actually  geometrical  but  is  included  here  because  of 
its  general  significance  and  because  the  tools  needed  to  prove  it  are 
available  in  this  appendix. 

PROPERTT  12;  Given  a  sequence  of  ellipses  to  be  combined  E^,  £2# 

if  none  were  rejected  then  the  sum  of  accept-reject  test  conditions  computed 
at  each  step 
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comes  out  the  same  regardless  of  the  order  of  the  E^.  (With  rejections 

the  result  does  not  have  the  form  shoun  above  because  of  missing  terms  and 
is  not  Independent  of  order.) 

Proof.  Proceed  by  by  induction.  For  n-2  the  symmetry  of  the  individual 
test  gives  the  result.  If  the  result  is  assumed  to  be  true  for  ngm  then 
the  n>m-^1  case  must  be  shorn.  Note  that  for  this  case  any  permutation  that 
does  not  move  the  m-»>1st  ellipse  would  have  the  same  k>m-^1st  term  in  the  sum  so 
that  it  is  sufficient  to  compare  the  value  of  the  sum  for  k-2  to  m.  But  since 
that  sum  applies  to  the  same  k  ellipses  our  Induction  hypothesis  applies. 

Thus  in  the  ease  where  the  ffl-»1st  ellipse  is  not  moved  our  induction  criteria 
is  satisfied. 

The  fact  that  the  induction  result  holds  when  the  position  of  the  last 
ellipse  is  not  permuted  will  allow  restriction  of  our  study  to  a  single  case. 
The  reason  for  this  simplification  follows.  If  the  result  holds  for  two 
permutations  a  and  B  then  it  holds  for  their  composition.  Thus  it  is 
sufficient  to  prove  that  the  generating  set  consisting  of  exchanges  of 
adjacent  ellipses  each  preserve  the  sum.  Only  one  of  these  permutations  moves 
the  last  ellipse.  Thus  it  is  sufficient  to  show  that  exchanging  the  last  two 
ellipses  preserves  the  sum. 

Then,  if  there  are  m  ellipses,  the  first  m-2  are  unchanged  when  the  last 
two  are  exchanged.  As  these  ffl-2  ellipses  are  constant,  they  may  be  combined 
into  a  single  ellipse.  The  resulting  expression  is  the  k-3  ease  with  the  new 
first  ellipse  being  the  same  as  the  old  combination  of  the  first  m'2  ellipses. 
Now  if  we  use  the  linear  transformation  sending 

to  the  origin  and  to  the  identity  matrix  then  it  remains  to  show  that 

x|(I+S2)"^X2+(09X2“X2)^(I®S2+Sj)“^ (OtXg-Xj)  is  unchanged  by  permuting 

ellipses  2  and  3.  First  compute  the  terms  in  this  expression. 

I®S2-(I+S“b“’-(S2^I)”^S2. 

0®X2-(  I®S2 )  (0+Sj^  X2 )- CSj*! 328"^  Xg-CSg+I  )"^  X2 . 

We  must  show  xT(l-.-S2)“^X2*(xJ(S2+I)"^*x|)((S2+I)"'‘S2+S2)’'’ ((S2+I)"‘‘X2*X2) 

-X^Cl+Sj)*^ Xg+CX^CSj^I)"^ -X|)(  (S3*I)"42^S2)"'' ( (Sj^I)"'’ X2-X2) 

Simplifying  this,  we  get  an  equation  of  the  form 

3^()X2+:^()X2+3^()X2*xJ()X2  -  XJ()X2^XJ()X2+X*()X2^X^()X2 

It  turns  out  that  the  corresponding  coefficients  in  this  quadratic  are  equal 
as  is  shown  below. 

Matching  x|()X2  type  terms  requires  that 

(I+S2)^!*(S2+I)')  ((S2^I)"'S2^S3)“^  (S2^I)“’-((S2*1)"'*S3+S2)"^ 

<->  (S2>I)"V(S2+S2(S2^I))“^(S2*I)"’-((S2*I)"^Sj+S2)“^ 

<->  (S2^S3(S2^I  ))♦!-(  (3 j^I )“’  ^  ^  ^  ^ 

<->  ( (Sj4l S2^S2) <S2^S2S2*S2*I )-(S2^I ) (S2+S3S2^S2) 

<->  (Sj^(S2*l)S2)(S2+S3S2+S3^I)-(Sj-*'l)(S2^I)CS2»S3S2+S3) 


When  Bultiplled  out  both  aides  have  12  terma  which  are  equal.  The  terms 
of  the  form  x|()X^  match  by  symmetry  with  the 
To  have  the  iJox,  terma  match  requires  that 

<->  ((S2^I)''S2^S2)')(S2+I)“'-(S3*I)”^(S3*I)"^S3*S2)'' 

<->  ((Sj*I)'''s2+S2)<S3*I)-<S2*I><<V^^’-V®3^ 

<->  S2+S2(S3^I)-S2^(S2*I)S3 

These  three  terms  match  when  multiplied  out.  The  terms  match  by 

symmetry.  Since  all  terma  match  the  expressions  are  equal. 


